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We discuss an explorative computation of real time autocorrelation functions, in 
the classical approximation. The results for the 'plasmon' frequencies and damping 
rates appear compatible with the divergencies expected from perturbation theory. 



1 Motivation 

The (effectively) classical approximation has been used in computations of the 
Chern-Simons diffusion rate in the SU(2)-Higgs model, which is an essential 
ingredient in theories of baryogenesis. It is potentially a very useful approxima- 
tion for non-equilibrium processes in field theory. However, its formal proper- 
ties are not very well known yet and it is not clear either how well the classical 
field theory approximates the quantum theory at high temperatures. 

These issues can be studied in the weak coupling expansion. The classical 
theory has divergencies at finite temperature and the question naturally arises 
if the theory is renormalizable in some sense. The answer is affirmative for 0"* 
theory and the hot classical theory can approximate the quantum theory well 
at weak coupling]. For gauge theories the situation is more complicated. Their 
leading hard thermal loop effects are non-local and counterterms would have 
to have this property as well (although the specific form of the nonlocalities 
may be representable by locally coupled new degrees of freedom). 

Still, one can suppose that any divergent quantity will be unobservable at 
low energies. Divergent frequencies and damping rates will be unobservable at 
sufficiently large time scales. The question is then, if there remain physically 
interesting phenomena at such time scales. 

For time independent quantities the classical approximation is a form of 
the dimensional reduction approximation. However, for the time component 
Aq of the^auge field the coefficients of Tr Aq and Tr Aq are zero in the classi- 
cal theorp, while they have well defined values in the dimensionally reduced 
theory, including a Debije mass counterterm cx TrAp. The absense of this 
counterterm in the classical theory implies that the corresponding Debije mass 
is divergent. For static quantities this is not a problem since Aq simply decou- 
ples. Such a decoupling approximation is often made in addition to dimensional 
resuction because the perturbative Debije mass is relatively large. For dynam- 
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ical properties one may hope that a similar decouphng is possible in the sense 
describe above. 

The situation is more serious in the dynamical case since also the classical 
plasmon frequency associated with the modes Ak, k = 1,2,3, is divergent. 
This can be seen in a one loop calculation on the lattice. The two phenomena 
(diverging Debije mass and plasmon frequency) are related by the fact that 
there is no local counterterm available for Af^ in the classical lagrangian. 

So it seems a good idea to compute nonperturbatively various correlation 
functions in the classical model and study their lattice spacing dependence. 
The results can then also be compared with analytic calculations, in the clas- 
sical and quantum theory. This work was undertaken with Wai Hung Tan^. 
We concentrated on autocorrelation functions 

Coit) = (0(00(0)) (1) 

of gauge invariant observables with the quantum numbers of the Higgs and W 
fields, at zero momentum, 

O^H = ip^ip, -^W^ ^iTT(bHdk-iA^itf3)cj)t^, (2) 

with if the Higgs doublet, (j) its matrix version {iT2ip*,ip) and ta ~ Tq/2 the 
SU(2) generators. In the Higgs phase these fields correspond to the usual Higgs 
and gauge fields. Since (/) — \/HV with V € SU{2), is essentially the 
gauge field in the unitary gauge = 1, as can be seen in perturbation theory 
by replacing if by its vacuum expectation value, neglecting fluctuations. In 
the plasma phase this reasoning breaks down as the theory behaves like QCD 
and perturbation theory runs into infrared problems. 



2 Simulation 

The effective clajSsical hamiltonian on a spatial lattice can be given schemati- 
cally in the fornH 



ff 



/^EtJ ^^KnA^ + ^^4^x + 5^ (1 - i Tr f/™„^) 



rj-} I- / ^ L 2 /nx 7nx ' " X ^ ' / J \ 12 

X ni<in 

+ (An<^x)^An'^x + A(<^i(^x - v'^ f]. (3) 

We use the parametrization and numerical implementation of KrasnitB. Com- 
parison with dimensional reduction leads to the following identification of 
parameters 



P ^ 4^^ (4) 



2 



= (,^8 

1 /149 + 9p 27 + 6p-p2 aT 27r] + 6pf] - p^fjV 
\ 486 ^ 27 1 27 jj ' 

where a is the lattice spacing, 77, ry and E are numerical constants, p = 
m'jj/m'^, and g^, mj^ are running couplings and masses in the 4D quantimr 
theory in the MS scheme with scale parameter p, = ATTe~'^T « 7T. Further- 
more, ze / = 1 to a, good approximation, which leaves only one parameter 
z ^ ze = Ztt to set the time scale relative to the lattice spacing. Intuitively 
one expects z = 1 + 0{g^, A), but recent work suggests that z may have to be 
taken substantially different from 1, depending on the lattice spacing. 

We used z — 1, g"^ — 4/9, p = 2A = 1, and various (3, v"^ and lattice sizes 
N^. First we did a scan with = 24 and fixed, implying varying T jmH by 
eq. (|). The ratio T/mn can be converted to T/T^ since the transition point 
Tc/iriH is known reasonably well. Then we spent more effort in the plasma 
phase with {N,f3) = (20, 11) and (32, 18.3), while was adjusted such that 
T/Tc = 1.52 in both cases. Hence, the lattice spacing changed by a factor of 
11/18.3 — 0.60 while the physics remained the same according to (||), including 
the physical volume: LT = NaT = 9N//3 ~ const. In the Higgs phase similar 
effort was spent with {N,/]) = (20,13) and (32, 21.7), without adjusting w^, 
however, which led to somewhat different T/Tc- 0.85 and 0.78, respectively. 

Fig. H shows examples of correlation functions. On the left we see Cnit) 
in the Higgs phase for /3 = 13 and 21.7. The damped oscillations can be fitted 
to the pole dominance form 

i?e~'''* cos(ijji + a) + const, (6) 

with a, const w 0. The plot for Cwit) looks similar, but with more damping 
and the additive constant effectively nonzero. In the plasma phase the results 
are qualitatively different. There Cnit) still shows oscillations with physical 
[auj < 1) frequencies, but with a much larger effective damping coefficient 7, 
and the additive background drops in time. For Cw{t) we find only a few 
oscillations at unphysically short times {t/a of order 2), after which a slow 
roughly exponential decay sets in (w — > and small 7 in (^)), shown on the 
right in Fig. |l|. Fig. ^ shows the results for lj and 7 in units of g^T (in lattice 
units, all the ato and 07 are reasonable small compared to 1). 

3 Interpretation 

Consider first the frequencies ujh,w in the Higgs phase, where H and W corre- 
spond to the usual Higgs and W fields in leading order perturbation theory. At 
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Figure 1: Left: Cuif) in the Higgs phase for /3 = 13 (small amplitude) and 21.7 (large 
amplitude). Right: plot of \nCw(f) in the plasma phase. The crossing lines correspond to 
/3 = 11, 18.3 at T/Ta = 1.52. The top curve corresponds to /3 = 11 at a different v'^ such 
that T/Tc = 1.23. The time is in lattice units: 'time' = t/a. The insert shows In Cu,(t) 
versus In t/a. 



first sight the smaU deviations of the 32'^ data from the other data (interpolated 
by eye) suggest lattice spacing independencep. This is what we optimistically 
highlighted in the first version of the paper □. However, our understanding 
increased substantially since then and as a result (being pressed as well by the 
referees of our paper) we reanalyzed the data, in the spirit of the expected 
^l\fa divergence in ujy/ from one loop perturbation theory. The upshot (ex- 
pressed in the revised version) is that the data are in fact compatible with 
this divergence. A recent analysis skows that also their magnitude can be 
understood from perturbation theory 

The damping rates on the other hand are expected to be finite to leading 
nontrivial order. In the Higgs phase they compare well with the analytic re- 
sults in the quantum theory. For example, the ^ = 21.7 {T/Tc = 0.78) data 
give -iwlg^T = 0.22(4) ^nd -ih/9^T = 0.028(4), to be compared with 0.176 
and ~ 0.02, respectivelyQ. In the plasma phase there are strong fluctations in 
Cnit) causing large errors on 7//; the errors on 7n/ are reasonable. In both 
phases the damping rate data are compatible with lattice spacing indepen- 
dence, which, given the large errors is not a very strong statement. However, 
the a independence of in the plasma phase {jw/g^T = 0.31(2), 0.028(4) 
for respectively (3 = 11, 18.3) suggests that it is a physical quantity as well. 
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Figure 2: 'Plasmon' frequencies uj/g'^T (left) and damping rates -i/g^T (right) versus T/Tc- 
Solid symbols correspond to H, open symbols to W\ circle: 24'^ data; diamond: 20'^ data; 
square: 32^ data. 



A more complete understanding may be expected from analysis of the full 
autocorpelation functions in perturbation theory. Work in this direction is in 
progress. For example, the qualitative. difference between Cnit) and Cw{t) at 
short times can be explained this wajO. A good understanding of the (over-?) 
damping in Cw{t) at larger times may also be obtainable this way. 
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